Theory of Fermi-surface pockets and correlation effects in underdoped YBa2Cu306.5 
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The detection of quantum oscillations in the electrical resistivity of YBa2Cu306.5 provides direct 
evidence for the existence of Fermi surface pockets in an underdoped cuprate. We present a theoret- 
ical study of the electronic structure of YBa2Cu307_5 (YBCO) aiming at establishing the nature of 
these Fermi pockets, i.e. Cu02 plane versus CuO chain or BaO. We argue that electron correlation 
effects, such as orbital-dependent band distortions and highly anisotropic self-energy corrections, 
must be taken into account in order to properly interpret the quantum oscillation experiments. 



The independent-particle picture fails miserably in 
the description of the undoped parent compound of the 
HTSCs, such as for instance La2Cu04 or YBCOg.o- 
These systems, as shown schematically in Fig.[T^ for a 
single square Cu02 plane or in Fig[2K,b for YBCOg.o, are 
predicted by band theory to be 1/2-filled metals; instead, 
due to strong electronic correlations, they are charge- 
transfer gap antiferromagnetic insulatorsi*^. While on 
the overdoped side of the phase diagram a well-defined 
band-structure-like Fermi surface is recovered, as in the 
case of heavily overdoped Tl2Ba2Cu06+^^'^'^, the under- 
doped regime is the most debated, with two alternative 
scenarios being proposed. The first originates from a 
wide variety of theoretical calculations^'^'^, which pre- 
dict that the first few holes doped in to the antiferro- 
magnet would give rise to four 'nodal Fermi pockets' 
with a volume, counting holes, proportional to the dop- 
ing X away from half filling (Fig.[T]3). The second con- 
sists of the truncation of the single-particle Fermi sur- 
face giving rise to four 'nodal Fermi arcs' (Fig.[T];); it is 
suggested based on experimental observations by angle- 
resolved photoemission spectroscopy (ARPES)^'^° and 
is naturally connected to the existence of a d-wave-like 
pseudogapiiii^ii^ii^. 

The discrimination between the two scenarios of 
Fig.[T|D,c has mainly relied on ARPES studies; one should 
realize, however, that this is a formidable and yet unset- 
tled taskl— >i^. It rests on the detection of an arc profile, 
or lack thereof, on the outer side of the pockets, where the 
electron-removal spectral weight is expected to be vanish- 
ingly small because of the strong drop in the quasiparticle 
coherence beyond the antiferromagnetic zone bound- 
ary (diamond in Fig.[TjD)ii. The discovery of a closed 
Fermi surface contour in quantum oscillation experiments 
on ortho-II YBCOe.s^^ might lead to a resolution of this 
longstanding conundrum, by providing direct support to 
the nodal Fermi pocket scenario of Fig.[T]D. For this as- 
signment to be conclusive, however, one must exclude 
the possibility that bands other than the Cu(3c?;j;2_j,2)- 
0{2px^y) from the Cu02 planes, such as for instance the 
BaO bands, are responsible for the detected signal. Note 
in this regard that the BiO and TIO pockets predicted 
for instance for B:^^ and TP" cuprates have never been 
detected experimentally; on the other hand, the situation 
for YBCO is much more complicated because of the pres- 



ence of the CuO chains, which can potentially give rise to 
Fermi surface sheets of purely chain character and pro- 
vide an additional channel for BaO hybridization. The 
presence of such BaO / CuO-chain based Fermi pockets in 
band theory calculations for optimally doped YBCO has 
been pointed out previously by Andersen et al.^^. 

To illustrate the band structure effects on the Fermi 
surface of the YBCO family, we have carried out den- 
sity functional theory calculations of the electronic struc- 
ture of YBCO7-5, for (5 = 0, 0.5, and 1 (Fig.[3 and 
Ref.[23). In agreement with previous density functional 
studie o^^'^^ , the Fermi surface of stoichiometric YBCOg.o 
and YBCO7.0 consists of two large cylindrical sheets, 
originating from the CUO2 bilayers, and small pockets of 
mainly BaO character (see below), all centered around 
the Brillouin zone corner (Fig.[2^-d). In addition, in 
YBCO7.0 there is also the open one-dimensional Fermi 
surface of the 'full' CuO chains running along the h- 
axis (Fig.[2jl). In oxygen ordered ortho-II YBCOe.s, the 
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FIG. 1: (Color online) Different scenarios for the low-energy 
electronic structure of the single, square CUO2 plane; Bril- 
louin zone dimensions are expressed modulo the lattice con- 
stant, (a) Independent-particle Fermi surface for a metal- 
lic, 1/2-filled Cu-0 band (green/gray); its volume is precisely 
50% of the zone. When correlations are considered, the sys- 
tem becomes an antiferromagnetic insulator; (b) light hole- 
doping might result in four Fermi pockets (red/gray) with 
a volume that, counting holes, corresponds to the doping x 
away from 1/2-filling (the enclosed diamond is the antiferro- 
magnetic Brillouin zone) . (c) Alternatively, the Fermi surface 
is truncated near the zone edges reducing to four disconnected 
arcs (blue/gray). In all panels, white and yellow (gray) ar- 
eas identify hole and electron momentum-space regions, re- 
spectively; the thickness of the Fermi surface lines represents 
the quasiparticle coherence Zi^ (i.e., the spectral weight of 
the electron-removal quasiparticle peaks), which is strongly 
fc-dependent in (b) and (c). 
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FIG. 2: (Color online) Local-density approximation (LDA) calculations of the in- plane electronic structure {k^ = 0), and 
corresponding two-dimensional projected Fermi surface, for (a,b) YBCOe.o, (c,d) YBCO7.0, and (e,f) oxygen-ordered ortho-II 
YBCOe.s (the width of the Fermi surface lines is proportional to the amount of dispersion), (g) Ortho-II YBCOe.s crystal 
structure, with full and empty CuO chains alternating along the a-axis and running parallel to the 6-axis (plane-octahedra 
and chain-plaquettes each contain one Cu atom), (h) Sketch of 1/4— Brillouin zone for ordinary YBCO and ortho-II YBCOe.s; 
because of the folding to the reduced zone, (7r/a,0)-(0,0) and (n /a,TT /b)-{0,-K /b) are sets of equivalent symmetry points in 
YBCOe.s. In (a-f), as throughout the paper, momentum-space coordinates are given modulo the lattice parameters a and b. 



unit cell is doubled along a-axis due to the alternation 
of full and empty CuO chains and, correspondingly, the 
Brillouin zone is folded along the (7r/2, 0) — (7r/2,7r) line 
(Fig.[2^,h; momentum-space coordinates are here given 
modulo the lattice constants). This affects both shape 
and number of Fermi surface sheets, but the small BaO 
pockets are still present although translated to (0,±7r). 

As one can see in Fig.[5f, if taken quantitatively the 
size of the (0,±7r) pockets is vanishingly small and is 
not nearly enough to explain the quantum oscillations 
in the electrical resistance observed by Doiron-Leyraud 
et alJ^. However, it is also evident that in such a multi- 
band case the Fermi surface is quite susceptible to very 
small changes of the relative position of the bands and/or 
to a minimal rigid shift of the Fermi level itself, which 
can be induced by deviations in oxygen concentration or 
the presence of disorder. For instance (Fig[3]), for Fermi 
level position AEf = 0, —13.6, and —19.1 meV the area of 
the (0,±7r) BaO pocket is ^fe = 0.3, 3.01, and 5.22 nm"^ 
respectively, to be compared to Ak = 5.1 nm^^ found in 
quantum oscillation experiments on YBCOe.g^- Also 
the BaO pocket average band mass, mb — 0.77, 1.68, and 
1.56 mo for the three values of AEp, is in qualitative 
agreement with the experimentally observed cyclotron 
mass m* = 1.9mo, where mg is the free electron mass^. 

Before elaborating on the deeper significance of these 
findings, we would like to discuss the anticipated BaO 
assignment of the small Fermi surface pockets. As shown 
in Fig. m for YBCOe.s, a strong peak in the density of 
states (DOS) at energies just below the Fermi level be- 



longs to Cu 3dxy and 3dyz orbitals which are, in fact, 
strongly mixed with the apex O 2p and Ba empty or- 
bitals in a pdir configuration^!. Thus, BaO-Cuchain is 
primarily the character of the flat bands (Fig. [2^) and 
corresponding small pockets (Fig.[2f) about (0,±7r) in 
YBCOg.s. A similar analysis for YBCOe.o confirms the 
same origin for the (±7r,±7r) pockets; as for YBCO7.0, 
which is somewhat of a different case because the CuO 
chains are all full, we find that oxygen from the BaO 
plane contributes the most to those bands defining the 
small hole pockets. Consistent with our assignment, all 
these pockets are rather isotropic in and ky, contrary 



AEF(meV) = AEF(meV) = -13.6 AEF(meV) = -19.1 






FIG. 3: (Color online) LDA Fermi surface of ortho-II 
YBCOe.s for different positions of the Fermi level, so as to 
show the rigid-band-like doping dependence of Fermi surface 
and, in particular, BaO-Cuc(iain pockets at (0,7r) and {tt,tt). 
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FIG. 4: (Color online) LDA density of states (DOS) for ortho- 
II YBCOe.s- The high DOS at Ep originates from BaO and 
(oxygen-empty) Cuchain orbitals; these are pdn orbitals with 
Cu dxz,yz character, thus very different from the pda orbitals 
typically found for the Cu02 planes (d^2_y2) or that one 
might expect for the empty CuO chains (d3z2_^2)^. 
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FIG. 5: (Color online) Density of states (DOS) and band dis- 
persion calculated for YBCOs.o within the LDA-I-U approxi- 
mation, with parameters U = 8eV and J^r = 1.34eV— The 
system is now a 1.08 eV gap charge-transfer insulator, with a 
Fermi level DOS dominated by the BaO-Cuchain bands; the 
CuO 2 plane band is split and pushed to high energies. 



to what is expected for a purely one-dimensional chain 
electronic structure. Also, they are very sensitive to the 
Cuchain-0 apex and -Ba distances, which show a substan- 
tial variation with S in YBCOy-a^. 

At this point the key question is: what is the relevance 
of the results we have just presented to the physics of 
the underdoped cuprates? After all, these are strongly 
correlated electron systems and it is well-known that the 
standard local-density approximation (LDA) to the ex- 
change correlation potential fails to treat strong corre- 
lations properly. In order to gain more direct insight 
into this very fundamental methodological shortcoming, 
we further our analysis within the so-called LDA+U 
method22i which, by adding the on-site Coulomb repul- 
sion U to LDA in a mean field-like way, can reproduce the 
correct magnetic ground state and forbidden gap in corre- 
lated antiferromagnetic insulators such as La2Cu04^'^^. 
However, the LDA+U approach works fairly well only in 
the case of integer orbital occupation and well-established 
long-range order (e.g., spin, charge, or orbital order)^^. 
This constrains our DFT analysis of the Fermi surface 
of YBCO in the presence of strong correlations to the 
single case of YBCOg.o (Fig-El) > which is truly an antifer- 
romagnetic insulator with 3cP and Sd^" configurations for 
plane-Cu^"'" and chain-Cu^"'' , respectively. It is important 
to realize that in such a framework, as shown in Fig.O 
the Cu02 in-plane states are almost completely removed 
from the Fermi energy and the first electron-removal state 
has nearly 100% BaO-Cuc/iam character. This is partly 
due to the inadequacy of LDA-I-U in fully accounting for 
the many-body physics of the Cu02 planes; while the 
LDA in-plane states are about 0.5 eV below the Fermi 
level, we do know that in the undoped cuprates there are 
Zhang- Rice singlet (ZRS) states"^° pushed ~0.5eV higher 
in energy, thus back into the first ionization energy range, 
but these two-particle states are completely missed by 



the single-particle LDA description^ . In the scope of this 
paper, we can make one very important observation: in- 
dependently of the exact energy of the ZRS states, in 
YBCOe.o a rigid-band-like hole doping would give rise 
also to pockets originating from the BaO-Cuchain bands 
and, in direct contrast to materials like La2Cu04, not 
just the in-plane Cu{3dj.2_y2)-0{2px^y) states. Further- 
more, the folding of the Brillouin zone, in this case due to 
antiferromagnetic spin order, translates the BaO-Cuc/iam 
pocket from {tt,tt) to the F-point (Fig.[T^ and[5)D). 

Let us now attempt a critical summary of experimental 
and theoretical findings. In light of our understanding of 
the underdoped cuprates, it is definitely very tempting 
to interpret the recent quantum oscillations results^^ as 
direct evidence for nodal Fermi pockets originating from 
the Cu02 planes electronic states. Our analysis shows, 
however, that the details of the band structure cannot 
be neglected in the case of YBCOe.s, since BaO-Guchain 
states might give rise to additional Fermi pockets con- 
sistent with the recent experimental observations both 
as far as volume and especially small effective mass are 
concerned. As a matter of fact, their role might become 
even more significant when correlation effects are also in- 
cluded, as they should for these very low doping values. 
Although not as extreme as in the case of YBCOg.o, the 
Coulomb repulsion U associated with the Cu 3d orbitals 
would strongly affect the in-plane Cu02 bands but much 
less so the BaO-Cucham bands. This would lead to an 
orbital-dependent band distortion (as opposed to a rigid 
shift) and trigger a hole transfer from the more to the 
less correlated electronic states, conspiring to generate 
qualitatively different Fermi features: 'highly-correlated' 
Cu02 pockets playing a primary role in the emergence 
of high-temperature superconductivity, and 'more con- 
ventional' BaO-Cuc/iain pockcts behaving as mere spec- 
tators. In the event of such coexistence, the latter would 
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be more likely observed in quantum oscillation experi- 
ments, while the Cu02 pockets might escape detection 
because of the highly momentum-dependent self-energy 
effects, which could even result in vanishing quasipar- 
ticle residues and/or diverging effective masses at some 
momenta along a Fermi pocket contour (e.g., beyond the 
antiferromagnetic zone boundary as in Fig. [Us). If this 
was indeed the case, it would also suggest that the de- 
tection of quantum oscillations might be limited to those 
cuprates containing CuO chains-^. 

Due to the complex multiband and correlated charac- 
ter of the electronic structure of YBCOe.a, the precise 
determination of the nature of the Fermi surface pock- 
ets revealed by quantum oscillation experiments^^ will 
require further study of this underdoped system by a va- 
riety of techniques. One serious issue, which has not been 
addressed here, concerns the negative sign of the Hall co- 
efficient observed by Doiron-Leyraud et al.^^. In a sim- 
ple one particle picture, this negative sign would point to 



an electron rather than hole-like character for the Fermi 
surface pockets detected in quantum oscillations. As in 
Ref.m, we also concentrated on the oscillatory part of 
the Hall resistivity above and beyond the rather smooth 
"background" contribution from which the negative sign 
is derived. We note that this Hall resistivity background 
is expected to be rather complicated because of contribu- 
tions from the strongly correlated Cu02 planes and the 
presence of localized magnetic moments; a net negative 
sign for the Hall coefficient need not be inconsistent with 
the presence of hole-like Fermi surface pockets. 
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Note: After completion of this work, we became aware 
of a recent DFT study of YBCO by Carrington and 
Yelland'^-. The LDA band structure results are very sim- 
ilar, with minor exceptions due to the difference in the 
treatment of the exchange correlation potential^^. 
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